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Procede de tatouage d'une image en couleurs par approche vectorielle, 
procede de detection d'une marque de tatouage, dispositifs, image et 
programmes d'ordinateurs correspondants. 

Le domaine de l'invention est celui du tatouage d'images fixes en 
5 couleurs. Plus precisement, l'invention concerne une technique de tatouage 
permettant d'integrer des signatures a une image, en tenant compte de la 
dimension couleur de cette derniere. 

Le domaine du tatouage d'images, encore appele « watermarking », est 
actuellement en plein essor, et fait l'objet de nombreuses recherches, tant dans le 
10 domaine des sequences video que des images fixes. 

On distingue actuellement deux grandes families de tatouage d'images, 
correspondant respectivement au tatouage d'images en niveaux de gris et 
d'images couleurs. 

En effet, les premieres annees de recherche dans le domaine du 
15 « watermarking » ont ete principalement axees sur des techniques de protection 
des images en niveaux de gris, et ont conduit a l'utilisation de trois domaines de 
marquage distincts : le domaine spatial, le domaine frequentiel et le domaine 
multiresolution. 

Les modeles de tatouage dans le domaine spatial (ou la marque est 
20 directement appliquee aux valeurs des pixels) ont l'avantage de necessiter un 
faible cout en temps de calcul. lis sont generalement robustes aux attaques 
geometriques (par exemple rotation et changement d'echelle). Des methodes 
proposees sont la modification d'histogramme (Coltuc D. et al., "Image 
authentication by exact histogram specification" (en francais "Authentication 
25 d' image par specification d'histogramme exacte"), workshop on multimedia 
signal processing, Cannes, France, octobre 2001) ou encore les techniques de 
patchwork (D. Gruhl, W. Bender, Moritomo, ''Techniques for data hiding" (en 
fran?ais, "technique de dissimulation de donnees"), in processing SPIE, volume 
2420, page 40, fevrier 1995). 
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Les modeles de tatouage dans le dornaine frequentiel ont quant a eux pour 
avantage d'etre robustes a la compression (de type JPEG par exemple). La marque 
est appliquee aux coefficients resultants d'une transformation du type Fourier 
(comme presente par exemple par V. Solachidis et I. Pitas, "Self-similar ring 
5 shaped watermark embedding in 2-D DFT domain 10 th European Signal 
Processing Conference EUSIPCO'2000, Tampere, Finlande, pages 1977-1980, 
septembre 2000) ou Cosinus Discret (comme presente par exemple par F. Alurki 
et R. Mersereau, "A robust digital watermark procedure for still images using 
DCT phase modulation 10 th European Signal Processing Conference 

10 EUSIPCO'2000, Tampere, Finlande, pages 1961-1964, septembre 2000). 

Le marquage dans le dornaine multiresolution offre enfin plusieurs 
avantages. II est avant tout le dornaine utilise dans les standards de compression 
les plus recents. II permet egalement de choisir la bande frequentielle qui porte la 
marque, permettant ainsi de reduire les risques de degradation de P image par 

15 Papplication de la marque (comme illustre par D. Kundur et D. Hatzinakos, 
t( Digital watermarking using multiresolution wavelet decomposition", 
Proceedings of IEEE ICASSP '98, vol. 5, pages 2969 - 2972, Seattle, WA, USA, 
Mai 1998). 

Outre ces techniques de tatouage des images en niveaux de gris, une 
20 deuxieme grande famille de techniques de tatouage d' images fixes propose de 
prendre en compte la dimension couleur des images. 

Au scin de cette famille, on distingue tout d'abord une premiere sous- 
famille de techniques consistant a adapter les methodes en niveaux de gris sur les 
trois composantes couleur. Des parametres sont alors utilises pour controler la 
25 force de marquage sur chaque composante afin de tenir compte des proprietes du 
systeme visuel humain. 

Une deuxieme sous-famille regroupe des methodes specifiques aux images 
couleur. Elle tient compte du systeme visuel humain et utilise les caracteristiques 
des representations de la couleur. 
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Par exemple, une methode particuliere, specifique aux images couleurs, est 
celle proposee par L. Akarun, N. Ozdilek, B. U. Oztekin, "A Novel Technique for 
Data Hiding in Color Paletted Images", Proceedings of the 10th European Signal 
Processing Conference, EUSIPCO'00, Tampere, Finlande, pages 123-126, 

5 septembre 2000. 

La premiere etape de cette methode consiste a quantifier l'espace de 
couleurs par l'algorithme dit de 'median-cut'. II consiste a partitionner de faeon 
iterative l'espace de couleur avec des plans perpendiculaires aux axes des couleurs 
et passant par les valeurs medianes des donnees. 

10 Deux cas sont alors possibles : 

- toute la palette est utilisee dans la representation de l'image, ou 

- certaines valeurs de la palette ne sont pas utilisees par l'image. 

Le premier cas n'est alors pas adapte a l'algorithme de marquage. Mais la 
palette obtenue comprend des couleurs que l'ceil humain ne peut discerner. En 
15 utilisant cette propriete on peut ainsi liberer certaines couleurs de la palette, de 
facon qu'elles ne soient pas utilisees dans l'image a marquee. 

L'auteur propose que deux couleurs sont indiscernables si AE < 3 avec 
AE = t/(L, - L 2 ) 2 + {a'-a 2 f +{b x -b 2 f , dans l'espace Lab. On rappelle que 
l'espace de representation Lab est un espace perceptuellement uniforme. L 
20 represente la luminance et les composantes a et b sont chromatiques. 

Soit M(i) la marque, composee de couleurs n'appartenant pas a la palette. 
L'auteur precise qu'une marque « binaire » (composee de deux couleurs) est plus 
robuste (le risque d'erreur de detection est ainsi diminue). 

Soit C la couleur la plus utilisee de la palette. (C est definie par le pic le 
25 plus eleve de l'histogramme couleur, chaque abscisse correspondant a une couleur 
de la palette.) 

La marque contient un nombre d'elements inferieur au nombre de pixels 
de couleur C, i<h(C). Chacun de ces elements est indiscemable de la couleur C. 
Le marquage consiste a remplacer le i 4me pixel de couleur C par M(i). 



1 er depot 

4 



Une autre methode originate est celle proposee par S. Battiato, D. 
Catalano, G. Gallo, R. Gennaro, dans "Robust Watermarking for Images based on 
Color Manipulation", Proceedings of the 3rd Workshop on Information hiding, 
LNCS 1768, pages 302-317, Dresden, 1999. Selon cette methode, la marque n'est 
pas creee pour l'image au prealable, mais c'est le contenu colorimetrique de 
1' image qui represente la marque. Un inconvenient de cette technique est done 
qu'elle est tres lourde en donnees a conserves 

L'espace de couleur propose par Fauteur respecte deux proprietes : 

l'espace doit etre perceptuellement uniforme (comme le sont les 
espaces Lab et Luv) de fa?on a pouvoir assimiler une mesure de 
distance euclidienne a une difference de couleur pour le systeme 
visuel humain ; 

le passage vers cet espace, note ZC,C 2 , doit etre rapide, simple et 

sans pertes d'informations. 
Le domaine des couleurs opposees est defini ici a partir de RGB, de la fa9on 
suivante : 

RGB — LC,C 2 : LC X C 2 ^RGB: 

L = R + G + B; f/?=(L + C 2 -C,)/3; 

, C, = 2B - R - G; " « G - (L - C 2 )/3 

C 2 =/?-2G + 5. [b = (C { +Z.)/3. 

D'apres les auteurs, cet espace se rapproche le plus de la representation des 
canaux chromatiques du systeme visuel humain. 

Le marquage est realise de la fafon suivante. La marque est un vecteur 
M(n), «=/,...,&,..., TV ou Af est le nombre de couleurs de l'image et£ une couleur de 
l'image. Soit (LC l C 2 ) k la couleur associee a l'indice k, representee par un 

vecteur dans l'espace LC { C 2 . On selectionne un rayon dans une sphere definie 
autour du point de cbordonnee {LC x C 2 ) k de fa^on aleatoire. Pour chaque pixel 
(x,y) correspondant a la couleur le vecteur correspondant au rayon (determine 
precedemment) est ajoute au vecteur de couleur initial pour obtenir le vecteur 
marque {LC { C 2 ) k % . 



1er depot 



5 



Chaque couleur est ainsi marquee par addition d'un meme vecteur. 
L'image marquee est reconstruite en rempla^ant les couleurs originates par les 
couleurs marquees, en respectant les coordonnees des pixels de l'image. 

La marque est alors constitute de Tensemble des couleurs de r image 
5 originale. 

La detection est realisee par comparaison de l'image marquee avec les 
marques generees sur 1 'ensemble des images traitees. Elle s'effectue de la fa<jon 
suivante. II faut d'abord considerer que l'image sur laquelle la detection est 
effectuee possede le meme nombre de couleurs que l'image originale. 
10 L'algorithme de detection compare l'image a chacune des marques enregistrees. 
Couleur par couleur, on va chercher la marque la plus proche (ou les vecteurs 
differences entre les deux pixels de meme coordonnee sont les plus semblables). 
La marque correspondant a l'image marquee est celle ayant le plus grand nombre 
de vecteurs de difference entre les couleurs de l'image marquee et de l'image 
15 originale en commun. 

On notera que les differentes publications relatives aux techniques de 
tatouage couleur sont generalement consacrees a l'integrite de l'aspect visuel et la 
robustesse de ces techniques contre les attaques classiques est peu developpee. 

Generalement, pour des methodes dont l'algorithme de base peut etre 
20 applique aux images en niveaux de gris, le marquage suit la technique suivante : 

- passage dans l'espace transformee (coefficients d'ondelettes, coefficients 
de cosinus discret, etc.) 
application de la formule suivante : 

25 bu l w represente la transformee de l'image (ou la composante) marquee, 

/' la transformee de l'image (ou la composante) originale, M la marque, a le 
facteur de controle d'intensite de marquage, et ou i et j represented les 
coordonnees du pixel traite. 

Les algorithmes bases sur les caracteristiques des composantes couleur 

30 sont plus adaptes aux caracteristiques du systeme visuel humain. Le premier 
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article de Kutter ML, Jordan F. et Bossen F, ("Digital Signature of Color Images 
using Amplitude Modulation", Processings of SPIE storage and retrieval for 
image and video databases, San Jose, USA, volume 3022, numero 5, pages 518- 
526, Fevrier 1997) sur le « watermarking » couleur proposait de travailler sur la 
5 composante bleue du systeme RGB, pour laquelle Toeil humain etait le moins 
sensible. 

Cette idee a ete reprise par A. Reed et B. Hannigan, dans "Adaptive Color 
Watermarking", Proceedings of SPIE, Electronic Imaging, volume 4675, janvier 
2002. Ces auteurs proposent de travailler sur la composante jaune du systeme 

10 CMY (« Cyan Magenta Yellow » pour « Cyan Magenta Jaune », qui est un espace 
colorimetrique), car ils considerent que Fceil humain est moins sensible aux 
variations de couleur sur l'axe jaune-bleu. 

De fa§on generate, la composante bleue permet de mieux dissimuler la"* 
marque (mais avec une robustesse moindre) et la composante verte permet de?/ 

15 mieux proteger la marque contre des attaques (mais avec des degradations plus • 
visibles de Timage). Cela s'explique par le fait que le systeme visuel humain est , 
plus sensible aux variations dans le vert qu'aux variations dans le bleu. Le»r 
compromis entre invisibilite et robustesse de la marque depend done de la nature^ 
colorimetrique de la composante. 

20 Enfin, une derniere technique, proposee par J.J. Chae, D. Mukherjee, et 

B.S. Manjunath, dans "Color Image Embedding using Multidimensional Lattice 
Structures", Proceeding of IEEE International Conference on Image Processing, 
Chicago, Illinois, volume 1 , pages 460-464, octobre 1998, repose sur une 
approche du type vectorielle. 

25 Selon cette technique, la premiere etape du marquage consiste a faire une 

decomposition en ondelettes de 1' image initiale e;t de la signature (qui peut aussi 
etre une image). Un seul niveau de decomposition est effectue. On obtient ainsi 
les coefficients en ondelettes de l'image originale, notes (C Y , C a , C v )(x, y) , dont 

chaque composante correspond a la composante couleur de Tespace YUV, et les 
30 coefficients en ondelettes de la marque. L' utilisation de l'espace YUV (espace 
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utilise en video : Y est la composante de luminance, U et V sont des composantes 
chromatiques) permet ainsi une adaptation directe de cette technique aux 
documents video. 

Les coefficients d' ondelettes resultant de la decomposition de la marque 
5 sont quantifies en j? niveaux. On obtient ainsi un vecteur M , 
(M yy M v , M v )(*, y) pour une marque en couleur et M (x, y) pour une marque en 
niveaux de gris, contenant ($.) elements ou 1 < i < /3 . 

L'integration de la marque peut s'ecrire sous la forme : 

(C Y , C u >Cy). t \x,y) = (C y? C u , C y ) i (x 9 y) + aMiS;) 

10 ou a est le facteur de controle de la force de marquage. 

La detection de la marque est ensuite effectuee de la fa?on suivante. Apres 
T application de la transformee en ondelettes sur V image marquee, les coefficients 
resultants sont quantifies en /? niveaux. 

Pour estimer le vecteur le plus proche de celui utilise en tant que marque, 
15 on recherche la valeur de quantification du coefficient le plus proche de ceux de 
Timage initiale (la marque est alors detectee, element par element). 

L'image correspondant a la marque detectee pourra ainsi etre reconstruite 
par transformee en ondelettes inverse. 

Toutes les techniques de tatouage d'images en niveaux de gris ou d'images 
20 couleurs decrites ci-dessus presentent de nombreux inconvenients. 

Ainsi, les techniques de « watermarking » des images couleur basees sur 
des algorithmes en niveaux de gris ne permettent pas la prise en compte de la 
dimension couleur des images. De ce fait, elles ne permettent pas d'optimiser le 
compromis entre invisibilite et robustesse de la marque. 
25 Les techniques de tatouage des images couleur sont quant a elles tres 

complexes et cherchent toutes a repondre a un objectif d'invisibilite de la marque. 
En d'autres termes, elles ne sont pas suffisamment robustes vis-a-vis des attaques 
classiques de type compression ou filtrage. 

La technique proposee par Chae et al. decrite ci-dessus, bien qu'utilisant 
30 une approche vectorielle, presente quant a elle pour inconvenient de ne pas etre 
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robuste a la compression. En outre, cette methode de marquage ne permet pas 
d'envisager la securisation de documents. Enfin, cette technique ne tient pas 
compte de la dimension couleur de l'image, puisqu'elle affecte la meme force de 
marquage aux trois composantes de couleur de l'image. 
5 Enfin, aucune technique de Tart anterieur utilisant une decomposition en 

ondelettes n'offre un marquage aveugle, qui permette une detection de la marque 
ne necessitant pas V image originale. Par ailleurs, aucune de ces techniques n'offre 
une robustesse contre la compression JPEG, le filtrage median et I'ajout de bruit. 

L' invention a notamment pour objectif de pallier ces inconvenients de Tart 
10 anterieur. 

Plus precisement, un objectif de l'invention est de fournir une technique de 
tatouage d'images couleur permettant d'integrer une ou plusieurs signatures a une 
image couleur. 

En d'autres termes, un objectif de r invention est de fournir une technique 
15 de securisation de documents par insertion invisible d'une marque. 

Un autre objectif de l'invention est de mettre en oeuvre une telle technique- 
qui soit robuste a la plupart des attaques. Notamment, un objectif de l'invention, 
est de fournir une telle technique qui soit robuste a la compression (notamment de 
type JPEG), au filtrage (notamment par un filtre median) et a I'ajout de bruit./ 
20 L'invention a encore pour objectif de mettre en oeuvre une telle technique 

qui permette une detection aveugle, c'est-a-dire ne necessitant pas V image 
originale, de la marque dans une image. 

Encore un objectif de l'invention est de fournir une telle technique 
permettant de reduire la visibilite de la marque par rapport aux techniques de Tart 
25 anterieur. 

Ces objectif s, ainsi que d'autres qui apparaTtront par la suite, sont atteints a 
l'aide d'un procede de tatouage d'une image en couleurs, presentant trois 
composantes. 

Selon l'invention, un tel procede comprend une etape d'insertion d'une 
30 marque de tatouage, en au moins un point de l'image, selon une regie d'insertion 
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prenant en compte la position relative de trois vecteurs de composante, pour 

chacune desdites composantes respectivement, associes audit au moins un point. 

Ainsi, l'invention repose sur une approche tout a fait nouvelle et inventive 

du tatouage d'images couleur. En effet, l'invention repose sur une approche 
5 vectorielle, permettant, par rapport aux techniques de l'art anterieur, de prendre en 

compte la dimension couleur de 1' image. 

La marque de tatouage n'est done plus, contrairement aux techniques de 

l'art anterieur, appliquee de facon identique sur les trois composantes couleur de 

l'image ; en d'autres termes, alors que, selon l'art anterieur, on procedait a un 
10 marquage independant des trois vecteurs de composante de l'image, selon 

l'invention, ces trois vecteurs sont desormais pris en compte « en synergie » pour 

I' insertion de la marque. 

La mise en oeuvre d'une approche vectorielle permet, en termes de 
robustesse, d'obtenir des reponses plus stables que les techniques de l'art anterieur 
15 dans le domaine des ondelettes, telles que celle de Kundur et al. citee ci-dessus, 
comme technique de tatouage dans le domaine des ondelettes. Une telle approche 
vectorielle permet egalement de minimiser les variations induites par la marque 
dans le domaine couleur, et ameliore done sensiblement l'invisibilite de la marque 
par rapport aux techniques de l'art anterieur. 
20 Avantageusement, on selectionne, pour chacun desdits points concernes, 

deux vecteurs en tant que vecteurs de reference et un vecteur a marquer pour 
porter ladite marque de tatouage. 

On notera que l'invention consiste done, contrairement a l'art anterieur, a 
marquer des vecteurs de composante (i.e. des vecteurs couleurs par exemple), et 
25 non de simples coefficients de l'image. 

Preferentiellement, ledit vecteur a marquer est le vecteur central, se 
trouvant entre lesdits vecteurs de reference. 

De maniere avantageuse, on calcule les distances deux a deux entre lesdits 
trois vecteurs, lesdits vecteurs de reference etant ceux les plus distants Tun de 
30 F autre. 
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De maniere preferential le, on determine une frontiere entre lesdits vecteurs 
de reference, definissant deux zones associees respectivement aux valeurs binaires 
« 0 » et « 1 ». 

On cree ainsi deux demi-espaces, associes respectivement a chacun desdits 
5 vecteurs de reference, et correspondant a une yaleur « 0 » ou « 1 » de la marque 
que Ton souhaite inserer. 

Avantageusement, ladite frontiere est definie par la bissectrice entre lesdits 
vecteurs de reference. 

Selon une caracteristique avantageuse de l'invention, le marquage dudit 
10 vecteur a marquer comprend un Replacement eventuel dudit vecteur a marquer 
dans l'une desdites zones, selon la valeur binaire de marquage a appliquer. 

Ainsi, si Ton souhaite inserer sur le vecteur a marquer une marque de 
valeur « 0 » et que le vecteur a marquer se trouve dans la zone ou le demi-espace 
associe a la valeur binaire « 1 », on deplace son extremite de fa?on qu'elle se 
15 trouve dans r autre demi-espace, ou zone, associes a la valeur « 0 ». 

Si Ton souhaite inserer sur le vecteur a marquer une marque de valeur 
« 0 » et que le vecteur a marquer se^ trouve deja dans la zone associee a la valeur 
binaire « 0 », on peut egalement deplacer son extremite de fa<jon a se rapprocher 
du vecteur de reference se situant dans ce demi-espace, de fa^on a augmenter la 
20 robustesse de la marque (tout en restant bien entendu dans le demi-espace associe 
a la valeur binaire « 0 »). 

Selon une autre caracteristique avantageuse de l'invention, ledit 
deplacement a une amplitude variable, fonction d'une force de marquage choisie 
en fonction d'au moins une caracteristique locale de ladite image. 
25 Ainsi, en fonction des caracteristiques colorimetriques ou des 

caracteristiques de texture locales de Timage, on peut choisir d'inserer une 
marque plus ou moins visible et plus ou moins robuste, en fonction de la force de 
marquage choisie. Une telle force de marquage peut prendre par exemple une 
valeur comprise entre 0 et 1, une force nulle correspondant a une marque peu 
30 robuste et invisible, une force egale a 1 correspondant a une marque robuste mais 
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visible: 

Selon une variante avantageuse de l'invention, un tel procede de tatouage 
comprend egalement une etape de transformation en ondelettes de chacune 
desdites composantes de V image, et lesdits trois vecteurs de composantes sont 
5 determines, pour chaque point d'au moins un niveau de decomposition de ladite 
transformation en ondelettes, pour chacune desdites composantes respectivement. 

L'invention combine done a 1'approche vectorielle une decomposition en 
ondelettes de 1 'image. L'utilisation d'une telle transformee en ondelettes permet 
avantageusement a la technique de l'invention d'etre plus robuste vis-a-vis 
10 d'attaques du type compression JPEG, filtrage median et ajout de bruit. 

On notera que par « point », on entend ici une position pixelique dans un 
niveau de decomposition de l'image, ainsi qu'illustre par la figure 10, qui 
represente une image echelle E 101, et les trois images de detail D N l R 102, D N 2,R 
103, D N 3 ' R 104 pour la composante rouge d'une decomposition en ondelettes dans 
15 un systeme RGB. Le « point » 105 est defini par une meme position pixelique 
dans chacune des trois images de detail referencees 102 a 104. 

Preferentiellement, on prevoit au moins deux conventions de marquage 
d'un vecteur. 

Selon une caracteristique avantageuse de l'invention, on choisit comme 
20 convention de marquage pour une image donnee celle limitant les risques de 
conflit lors de la detection de ladite marque. 

Preferentiellement, on choisit une desdites conventions en fonction des 
nombres de vecteurs de reference present dans ladite image pour chacune desdites 
composantes. 

25 Ainsi, dans un systeme de type RGB par exemple, on determine le nombre 

de fois ou les vecteurs rouge, vert et bleu sont choisis comme vecteurs de 
reference ou comme vecteur a marquer, et on en deduit la convention de 
marquage permettant de reduire au maximum les risques de conflit lors de la 
detection de la marque. Par exemple, on choisit une convention de marquage telle 

30 que le vecteur le moins souvent marque soit situe dans la zone de conflit. 
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De maniere preferentielle, ledit marquage est, au moins dans des situations 
normales, calcule selon l'equation suivante : 

V„ jr ix,y) = V H (x.y) - (1 - F m )(V R (x 9 y) -V M (x,y)) 
oil V M (x 9 y) est ledit vecteur a marquer, V R (x,y) est un desdits vecteurs de 
5 reference, F M est ladite force de marquage et V M w (x 9 y) est ledit vecteur marque. 

Preferentiellement, en cas de conflit, ledit marquage est calcule selon 
l'equation suivante : 

V M , v (x>y) = V H (x,y)-beta.(l - F M )(V H {x,y) - V u {x 9 y)h 

ou beta<L 

10 Avantageusement, un tel procede de tatouage comprend, apres ladite etape 

d'association d'une marque de tatouage, une etape de transformation en ondelettes 
inverse, delivrant une image marquee. 

Ainsi, apres avoir procede a une transformation en ondelettes des ; 
differentes composantes couleur de 1' image, et apres avoir insere la marque de v . 
15 tatouage sur les vecteurs couleurs, on reconstruit une image marquee, par 
transformation en ondelettes inverse. 

De fa§on avantageuse, ladite marque est une signature binaire pseudo- 
aleatoire ecrite de fagon redondante. 

La redondance de Finsertion de la marque permet une meilleure protection 
20 de 1' image et une plus forte probability de detection du tatouage. 

Preferentiellement, lesdites composantes appartiennent au groupe 
comprenant : 

les composantes RGB ; 
les composantes YUV ; 
25 - les composantes CMY. 

Plus generalement, tous les espaces de representation de la couleur 
peuvent etre utilises, qu'il s'agisse d'espaces definis par la CIE (Commission 
Eclairage Internationale) ou d'espaces originaux. 

L'invention concerne aussi un dispositif de tatouage d'une image en 
30 couleurs, presentant trois composantes, et comprenant des moyens d'insertion 
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d'une marque de tatouage, en au moins un point de l'image, selon une regie 
d' insertion prenant en compte la position relative de trois vecteurs de composante, 
pour chacune desdites composantes respectivement, associes audit au moins un 
point. 

5 L'invention concerne egalement un programme d'ordinateur comprenant 

des instructions de code de programme enregistre sur un support utilisable dans un 
ordinateur pour le tatouage d'une image en couleurs, presentant trois 
composantes. Selon Finvention, ledit programme comprend des moyens de 
programmation lisibles par ordinateur pour effectuer une etape d'insertion d'une 

10 marque de tatouage, en au moins un point de l'image, selon une regie d'insertion 
prenant en compte la position relative de trois vecteurs de composante, pour 
chacune desdites composantes respectivement, associes audit au moins un point. 

L'invention concerne encore une image en couleurs presentant trois 
composantes et comprenant un tatouage obtenu par insertion d'une marque de 

15 tatouage, en au moins un point de l'image, selon une regie d'insertion prenant en 
compte la position relative de trois vecteurs de composante, pour chacune desdites 
composantes respectivement, associes audit au moins un point. 

L'invention concerne aussi un procede de detection d'un tatouage dans une 
image tatouee, realise selon le procede de tatouage decrit precedemment Selon 

20 l'invention, un tel procede de detection comprend une etape de recuperation d'une 
marque de tatouage, en au moins un point de l'image, selon une regie de 
recuperation prenant en compte la position relative de trois vecteurs de 
composante, pour chacune desdites composantes respectivement, associes audit au 
moins un point. 

25 Preferentiellement, ladite etape de recuperation comprend une sous-etape 

de calcul des distances deux a deux entre lesdits trois vecteurs, et les deux 
vecteurs les plus distants l'un de 1'autre sont des vecteurs de reference, le 
troisieme vecteur etant le vecteur marque portant ladite marque de tatouage. 

Avantageusement, une frontiere entre lesdits vecteurs de reference 

30 definissant deux zones associees respectivement aux valeurs binaires « 0 » et 
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« 1 », ladite etape de recuperation comprend egalement une sous-etape de 
reperage de la zone dans laquelle se situe ledit vecteur marque, et une etape 
dissociation d'une valeur binaire correspondante. 

De facpon avantageuse, ladite marque est recuperee au moins deux fois, et 
5 on met en oeuvre un calcul de correlation par rapport a une signature de reference, 
afin de decider si le tatouage est ou non correctement detecte. 

L'invention concerne encore un dispositif de detection d'un tatouage dans 
une image tatouee, realise selon le procede de tatouage decrit precedemment. 
Selon l'invention un tel dispositif comprend des moyens de recuperation d'une 
10 marque de tatouage, en au moins un point de Timage, selon une regie de 
recuperation prenant en compte la position relative de trois vecteurs de 
composante, pour chacune desdites composantes respectivement, associes audit au 
moins un point. 

L'invention concerne enfin un programme d'ordinateur comprenant des 
15 instructions de code de programme enregistre sur un support utilisable dans un 
ordinateur pour la detection d'un tatouage dans une image tatouee, realise selon le 
procede decrit precedemment. Un tel programme comprend des moyens de 
programmation lisibles par ordinateur pour effectuer une etape de recuperation 
d'une marque de tatouage, en au moins un point de Fimage, selon une regie de 
20 recuperation prenant en compte la position relative de trois vecteurs de 
composante, pour chacune desdites composantes respectivement, associes audit au 
moins un point. 

D'autres caracteristiques et avantages de l'invention apparaltront plus 
clairement a la lecture de la description suivante d'un mode de realisation 
25 preferentiel, donne a titre de simple exemple illustratif et non limitatif, et des 
dessins annexes, parmi lesquels : 

la figure 1 presente un synoptique du principe d'analyse multi-resolution 

d'une image I par transformation ondelettes, mis en oeuvre selon 

F invention ; 

30 - la figure 2 illustre le principe d'obtention d' images au niveau de resolution 
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j a partir de Fimage echelle au niveau j+1 par transformation ondelettes, 
selon le principe de la figure 1 ; 

la figure 3 decrit une arborescence de coefficients ondelettes construite a 
Tissue de la transformation ondelettes de la figure 2 ; 
' 5 la figure 4 illustre la generation d'une marque par redondance de la 

signature ; 

la figure 5 presente Fetape de construction des vecteurs de reference et du 
vecteur de marquage ; 

les figures 6a a 6c illustrent les differentes conventions de marquage mises 
10 en oeuvre selon Finvention ; 

• - ' la figure 7 illustre un cas de figure possible de positionnement des vecteurs 
originaux, ainsi que les deux deplacements du vecteur a marquer 
envisageables pour ce cas ; 

les figures 8a et 8b presentent respectivement un exemple d'image 
15 originate et d'image correspondante marquee selon la technique de 

tatouage de Finvention ; 

les figures 9a a 9c presentent, sous forme de graphes, les resultats de tests 
de robustesse de la technique de tatouage de Finvention, sur Fexemple de 
F image de la figure 8b ; 
20 - la figure 10 precise la notion de « point » d'une image, qui represente, 
dans le cadre d'une transformation en ondelettes, une position pixelique 
dans un niveau de decomposition. 

Le principe general de Finvention repose sur la prise en compte de la 
dimension couleur des images, et propose une approche vectorielle du tatouage. 
25 L' invention permet ainsi un marquage vectoriel prenant en compte, en synergie, 
les trois vecteurs couleur de Fimage. 

Dans un mode de realisation particulier de Finvention, qui sera decrit dans 
toute la suite du document, on combine a cette approche vectorielle une 
decomposition en ondelettes des composantes couleur de Fimage. 
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Dans ce mode de realisation particulier, 1' invention consiste done a inserer 
une marque dans une image couleur en realisant les etapes suivantes : 

- Transformation de chaque composante de 1' image couleur dans le domaine 
des ondelettes, 

5 - Generation independante de la marque de fa?on pseudo-aleatoire a partir 

d'une cle, 

- Definition de Tensemble des vecteurs couleurs 

Pour chaque coordonnee des composantes passe-bande: 

o Definition de 2 vecteurs de reference et d'un vecteur qui portera la 
10 marque (les vecteurs de reference sont les vecteurs les plus 

distants, la marque est ainsi appliquee sur le vecteur « central » ce 
qui minimise Timpact visuel) 
o Definition d'une regie qui gere la force de la marque et determine 
l'espace d'insertion associe au bit 1 et au bit 0. .t 
15 • - Reconstruction de l'image couleur marquee par transformee en ondelettes, 
en tenant compte des modifications par marquage des coefficients - 
. d'ondelettes. 

Ces differentes operations sont detaillees' plus en detail dans la suite, en f : 
relation avec les figures. 

20 1. Transformation ondelette 

La premiere etape mise fen ceuvre selon 1' invention est une etape de 
transformation en ondelettes. On rappelle que la transformation ondelette est un 
outil mathematique puissant permettant Tanalyse multi-resolution d'une fonction, 
ainsi que decrit par Mallat S. (dans "A Theory for Multiresolution Signal 

25 Decomposition: the Wavelet Representation", IEEE Transactions on Pattern 
Analysis and Machine Intelligence, Vol. 11, No. 7, Juillet 1989, pp. 674-693), par 
Stollnitz E.J., DeRose T.D., et Salesin D. (dans '"Wavelets for Computer 
Graphics: A Primer-Part V\ IEEE Computer Graphics and Applications, Mai 
1995, pp.76-84 ou encore dans "Wavelets for Computer Graphics: A Primer-Part 

30 2", IEEE Computer Graphics and Applications, Juillet 1995, pp.75-85). 
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L'annexe 1, qui fait partie integrante de la presente demande de brevet, 
rappelle les aspects essentiels de la theorie des ondelettes. 

Dans le cadre de 1' invention, les fonctions considerees sont des images 
numeriques, c'est-a-dire des fonctions bi-dimensionnelles discretes. Sans perte de 
5 generalites, on suppose que les images traitees sont echantillonnees sur une grille 
discrete de n lignes et m colonnes et a valeur dans un espace de luminance 
echantillonne a 256 valeurs. De plus, on suppose que n = 2 k (kGZ) et que 
m = 2 l (IfEZ). 

Si on note I Timage originate, on a alors : 

[0,wM0,»]-* [0,255] 
(x,y)v-*I(x,y) 

La transformation ondelette de r image I 10 permet une representation 
multi-resolution de I, ainsi qu'illustre par la figure 1. A chaque niveau de 
resolution 2 J (j<<>-l), la representation de I 10 est donnee par une image 
grossiere A^J 11 et par trois images de details D\,I 12, D\j 13 et D\j 14. 

15 Chacune de ces images est de taille 2 k * J x 2 UJ . 

A nouveau, au niveau de resolution suivant, Timage grossiere A y I 11 

peut etre representee par une image grossiere A ,_,/ 11, et par trois images de 

details D\^I 12,, D\ h J 13, et D^.J 14,. 

La transformation ondelette necessite le choix d'une fonction echelle 
20 ®(x) ainsi que le choix d'une fonction ondelette *V(x). A partir de ces deux 

fonctions, on derive un filtre echelle H et un filtre ondelette G dont les reponses 

impulsionnelles respectives h et g sont definies par : 

h(n) = ^ 2 _! ( a ),(p( u-n fpJnEiZ 

g( n) - ty 2 +( u)M u - «j)Vn£Z. 

Notons respectivement //et G les filtres miroirs de H et G (i.e. h(n) = h(-n) et 
25 g(n) = g(-n)). 

On peut alors montrer, comme illustre par la figure 2, que : 
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❖ A lf I 20 peut etre calculee en convoluant A^I 21 avec i/dans les deux 

dimensions et en sous-echantillonnant d'un facteur deux (23, 24) dans les 

deux dimensions ; 

❖ D\jl 22 peut etre calculee en : 

1. convoluant A^I 21 avec //suivant la direction y et en sous- 
echantillonnant d'un facteur deux (23) suivant cette meme 
direction ; 

2. convoluant le resultat de l'etape 1) avec G suivant la direction x et 
en sous-echantillonnant d'un facteur deux (25) suivant cette meme 

direction. 

❖ D*jl 26 peut etre calculee en : 

1. convoluant A 2 j+J 21 avec G suivant la direction y et en sous- 
echantillonnant d'un facteur deux (27) suivant cette meme 
direction ; 

2. convoluant le resultat de l'etape 1) avec //suivant la direction x et 
en sous-echantillonnant d'un facteur deux (28) suivant cette meme 

direction. 

❖ D^ I 30 peut etre calculee en : 

1. convoluant A jh J 21 avec G suivant la direction y et en sous- 
echantillonnant d'un facteur deux (27) suivant cette meme 
direction ; 

2. convoluant le resultat de l'etape 1) avec G suivant la direction x et 
en sous-echantillonnant d'un facteur deux (29) suivant cette meme 
direction. 

Dans un mode de realisation particulier de 1' invention, on utilise la base de 
Daubechies, et on choisit le niveau de resolution 2 r (r <s -\) avec r=-4. 

Les ondelettes de Daubechies ont un support de taille minimale pour un 
nombre donne de moments nuls. Les ondelettes de Daubechies a support compact 
se calculent a partir de filtres miroirs conjugues h de reponse impulsionnelle finie. 
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De la longueur du filtre depend le support temporel de Pondelette associee ainsi 
que le nombre de moments mils. II n'existe pas de forrnule explicite pour les 
filtres de Daubechies quel que soit l'ordre du filtre. Les filtres se construisent a 
partir de la resolution du probleme suivant : 

Soit le polynome trigonometrique h(w) = ^h[n]e~ ittw transformee de Fourier du 
filtre de Daubechies, alors il faut construire un polynome R(e~ tw ) de degre 
minimal tel que h(w) = yfli R(e~ iw ) et h{w) 



2 + 



h{w + 7l) 



= 2. Selon 



2 

T invention, on choisit de preference le filtre a huit coefficients qui correspond au 

polynome suivant : 

h{w) = 0,2304 +0,7 148^" iV +0,6309.e" <2vv -0,028.£-' 3 " -0,187^' 4vv 

+ 0,0308.<T 5m ' + 0,0329.<T r6,v -0,0106^ ,7w 

La figure 3 illustre la construction, apres que la transformation ondelettes 

est effectuee jusqu'a un niveau de resolution minimal 2 r (r<-l), d'une 

arborescence des coefficients ondelettes : 

- Chaque pixel p(x,y) reference 34 de Pimage A 2 ,I est la racine d'un arbre ; 

15 - Chaque racine p(x,y) se voit attribuer trois nceuds fils references 31 a 33 

designes par les coefficients ondelettes des trois images de details 
(s= 1,2, 3) localises au meme emplacement (x,y) ; 

- Du fait du sous-echantillonnage d'un facteur deux effectue par la 

transformation ondelette a chaque changement de resolution (voir figure 
20 2), chaque coefficient ondelette D* r (x,y) (s= 1,2,3) correspond a une zone 

de taille 2x2 pixels dans Pimage detail correspondante a la resolution 
2 r ^ . Cette zone est localisee en (2x,2y ) et tous les coefficients ondelettes 
y appartenant deviennent les noeuds fils de D s r (x y y). 

De fa9on recursive, on construit ainsi 1' arborescence dans laquelle chaque 
25 coefficient ondelette D s r (x,y) (s=l,2,3 et 0>u>r) possede quatre nceuds fils 

designes par les coefficients ondelettes de Pimage D^J localises dans la region 
situeeen (2x,2y) et de taille 2x2 pixels. 
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Une fois l'arborescence construite, chaque coefficient ondelette 
D s ,(x,y)(s=\,2, 3) correspond a une region de taille 2" *2~ r pixels dans Limage 



detail D[J 



2. Construction de la marque 

5 Une signature binaire S de N_ bits est generee de facon pseudo-aleatoire, 

controlee par une clef K. Cette signature, ecrite sous forme d'une matrice 
■y[N*-}fNest ensuite ecrite de facon redondante pour obtenir la marque W. La 
redondance peut-etre effectue bit par bit ou signature par signature, ainsi 
qu'illustre par la figure 4. 

10 Dans un mode de realisation particulier de 1'invention, la signature S 40 est 

une signature binaire de 16 bits et est generee de facon pseudo-aleatoire a 1'aide 
d'une clef K=9. Cette signature matricielle (4*4) 40 est ensuite recopiee T fois 
pour creer une marque, comme dans le premier exemple 41 presentee dans la 
figure 4. Pour une image de taille 256*256, sachant que les coefficients 

15 d'ondelettes utilisees sont ceux definis pour la 4 4me echelle de decomposition, la 
marque sera de taille 16* 16. La signature est done repetee 16 fois. 

On notera que sur la figure 4, les carres noirs illustrent un bit a 1 et les 
carres blancs represented un bit a 0. r 

3. Insertion de la marque 

20 3.1 Definition des vecteurs 

Les vecteurs sont definis a partir de la decomposition en ondelettes a 
Techelle -r, pour chaque coordonnees (x,y) des trois images de details D\jl 22, 
D 2 2 ,I 26, D\J 30, et pour chacune des trois composantes couleur notees I={A, B, 

C}. 

25 A, B et C sont des composantes qui dependent de l'espace de 

representation (ABC peut etre RGB, XYZ, YUV, etc.). 

On obtient ainsi, pour chaque coordonnees (x,y) 50 des images de details, 

les trois vecteurs suivants, illustres par la figure 5 : 

V A (x, y) = {D \, A(x, y),Djj A(x, y), D \, A(x, y) } , 
30 V B (x,y) = {D^B(x,y),D^B(x,y),D 3 2i B(x,y)}, 
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V c (x, y) = (D^C(x, y), D|,C(x, y), D^C(x, y) } 

Dans un mode de realisation particulier de P invention, on considere 

Timage dans l'espace de representation RGB, et les vecteurs sont definis pour 

chaque composante (R,G et B) a partir des coefficients d'ondelettes de la 4* me 

5 echelle de decomposition (r-4): 

V R (x,y) = {D^R(x.y), D^R(x,y), D^R(x,y) }> 
V G (x, y) = {D;,G(x, y),D*.<G(x, y), D 3 2 .G(x, y) } , 
V 8 (x, y) = {D 1 ^ B(x, y), D>., B(x, y), B(x, y) } 

3.2 Calcul des distances, vecteurs de reference et vecteur marque 
10 Pour chaque coordonnee des images details, les distances entre vecteurs 

deux a deux sont calculees : 

D A . B (x,y)-|V A (x > y)-V B (x,y)| 

D A;C (x > y) = |v A (x,y)-V c (x,y)| 

DB, c(^y)=|v B (x J y)-v c (x,y)| 

15 La distance la plus grande est definie par les deux vecteurs de reference 

V R 51 et V R ^ 52 (voir figure 5). Ce sera le troisieme vecteur V M 53 qui sera 

marque. 

Ainsi, si D AB (x,y)>D AC (x,y) et D AB (x,y) > D BC (x,y), alors 
V A (x,y)et V B (x,y) sont les vecteurs de reference et V c (x,y) est le vecteur qui 

20 portera la marque. 

Si D AC (x,y)>D AB (x,y) et D AC (x,y) > D BC (x,y), alors V A (x,y)et 

V c (x,y) sont les vecteurs de reference et V B (x,y) est le vecteur qui portera la 
marque. 

Si D BC (x,y)>D AB (x,y) et D BC (x,y) > D AC (x,y), alors V B (x,y)et 
25 V c (x,y) sont les vecteurs de reference et V A (x,y) est le vecteur qui portera la 
marque. 

Dans un mode de realisation particulier de l'invention, les distances entre 

les vecteurs sont calculees, deux a deux, pour chacune des 256 coordonnees 

possibles (16*16 coefficients d'ondelettes dans les images de details). 
30 . D R ^(x,y) = |V R (x,y)-V G (x,y)| 

■D RlB (x,y) = |v R (x >y )-V B (x,y)| 
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D G .B( x >y) = |Vc(^y)-V B (x,y)| 
Si D RiG (x,y)>D R>B (x,y) et D R0 (x,y) > D GB (x,y) , alors V R (x >y )zt 
V G (x,y) sont les vecteurs de reference et V B (x,y) le vecteur qui portera la 

marque. 

5 Si D R B (x,y)>D RG (x,y) et D R B (x,y) > D G B (x,y), alors V R (x,y)et 

V B (x,y) sont les vecteurs de reference et V c (x,y) le vecteur qui portera la 
marque. 

Si D G B (x,y)>D RG (x,y) et D G B (x,y) > D R B (x,y), alors V c (x,y)tl 
V B (x,y) sont les vecteurs de reference et V R (x,y) le vecteur qui portera la 
10 marque. 

3.3 Definition et choix de la convention de marqnage 
Pour chaque coordonnee des images de details, on cherche a connaitre le 
nombre de fois que les vecteurs V A (x,y), V B (x,y) et V c (x,y) sont utilises 

comme reference. Cette etape permet ainsi de choisir la meilleure convention pour 

15 1'operation d'insertion (voir figures 6a a 6c) 5 pour minirniser le conflit qui sera 

explique plus en detail au §3.5. 

Si V A (x,y) est le vecteur le moins utilise comme reference, on choisit la 

convention 3, illustree en figure 6c. 

Si V B (x,y) est le vecteur le moins utilise comme reference, on choisit la 

20 convention 1 , illustree en figure 6a. 

Si V c (x,y) est le vecteur le moins utilise comme reference, on choisit la 

convention 2, illustree en figure 6b. 

Ainsi, dans le cas de 1'image 'Maison' illustree en figures 8a (image 
originale) et 8b (image marquee), le vecteur V R (x,y) est utilise 235 fois en tant 
25 que vecteur reference, le vecteur V G (x,y), 108 fois et le vecteur V B (x 9 y), 169 

fois. Pour minirniser le conflit, on utilise done la convention 1, presentee sur la 
figure 6a. 

Toutefois, il est conseille de faire un choix precis dans le cas des images 
de petite taille. Le choix automatique de la convention de marquage est mieux 
30 adapte aux images de grande taille (au moins quelques milliers de pixels). 
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3.4 Operation d'insertion 

Le marquage consiste a modifier le vecteur en placant son extremite plus 
pres de l'extremite de l'un des deux vecteurs de reference, en fonction de la valeur 
de la marque (0 ou 1), et cela pour chaque coordonnees 50, ainsi qu'illustre par la 
5 figure 5. On note V RX {x,y) et V R2 (x,y) (51 et 52) les deux vecteurs de 
' reference et V M (x, y) le vecteur a marquer 53 . 

De facon generate, le marquage se presente sous la forme : 

avec suivant le cas : V R = V RX ou V R = V R2 suivant la valeur de M et la convention 
10 choisie. 

F M represente la force de marquage appliquee au vecteur V M (x,y) avec M={A, 
B, C}. La force de marquage peut etre ainsi adaptee en fonction des 
caracteristiques (colorimetrique ou de texture) locales de l'image. 

Deux cas de figure peuvent se presenter, ainsi qu'illustre pour un cas 

15 particulier par la figure 7 . 

Suivant la convention choisie et suivant la valeur de la marque W(x,y) : _ 

le vecteur V M (x, y) 73 est deja dans le demi-espace correct ou V M W (x, y) 

doit etre place. Dans ce cas, la formule (1) peut etre appliquee directement 
*V M (x,y)13. 

20 Ce cas de figure est illustre sur la figure 7 , en considerant le cas particulier 

ou W(x,y)=l, et ou la convention impose que V M }V (x, y) soit dans le demi- 

espace contenant V R2 (x,y) 72. 

le vecteur V M (x,y) 73 n'est pas dans le demi-espace correct ou 
f u w ( x ,y) doit etre place. Dans ce cas, avant d'appliquer la formule (1) a 

25 V M (x, y) 73 , il faut modifier la position de 1 'extremite de ce vecteur. Pour 
cela, il faut deplacer cette extremite dans une direction colineaire au 
vecteur V col (x,y) =V Rl (x,y)-V R2 (x,y) , d'une distance (D m +e), ou D m est 
equivalente a la distance entre V M (x,y)13 et (P M ) (ou P M est la frontiere 
entre les deux demi-espaces associes a chacun des vecteurs de reference71 

30 et 72). On choisit de preference e egal a 10% de D m , de facon a ce que 

V M ( x ,y) 73 ne soif pas confondu avec (P M ). 
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Ce cas de figure est illustre sur la figure 7, en considerant le cas particulier 
oil W(x,y)=0 9 et ou la convention impose que V M W (x y y) soit dans le demi- 

espace contenant V RX (x 9 y) 71. 

On notera que dans le cas ou F M = 0 , le marquage est de robustesse 
5 minimum. Le vecteur V M (x 7 y) 73 devient V M W (x, y)par simple deplacement (si 

necessaire) dans le demi-espace correspondant aux positions des vecteurs de 
reference 71 et 72, a la valeur de la marque et a la convention choisie. L' image est 

ainsi degradee au minimum. 

Dans le cas ou F M = 1 , la force de marquage est maximum. La marque est 
10 tres robuste, mais a une grande probability d'etre visible. V M W (x,y) devient alors 
confondu a Tun des vecteurs de reference 71 et 72 en fonction de la situation 

precedemment indiquee. 

Le fait de choisir une valeur intermediaire a F M permet de controler le v 

compromis entre invisibilite et robustesse de la marque. 
15 3.5 Conflit 

Les schemas de convention de marquage definissent le sens du marquage 
pour chacun des cas de figures possibles. Dans le cas ou la force de marquage ~ 
F„=l c'est-a-dire qu'elle est maximum, le vecteur marque devient egal a unr. 

At ' x 

vecteur de reference. 

20 Ainsi, sur les figures 6a a 6c , on remarque que, pour chacune des 

conventions, les bits de la marque sont opposes pour l'une des branches des 
triangles de convention. Un conflit est alors possible. En effet, lors de la 
transmission de Timage marquee ou des traitements qu'elle peut subir, les deux 
vecteurs confondus peuvent etre sensiblement modifies. II est alors possible que le 

25 vecteur de reference devienne le vecteur marque et vice- versa. 

Lors de la detection de la marque (qui sera decrite plus en detail au §5 ci- 
dessous), des erreurs de detection peuvent alors etre integrees a cause de ce 
conflit. Le conflit est minimise si la convention choisie minimise ce cas de figure. 
L'etape de choix de la convention proposee au §4.3 ci-dessus est done 

30 particulierement importante. 
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Dans le cas d'un conflit, 1'operation de marquage devient : 
V u w (x,y) = V K (x,y) - beta.Q. - F M ){V K (x,y) - V M (x,y)) 
oil beta<l. 

Ainsi le choix de la convention permet de minimiser le nombre de conflits, 
5 et les erreurs dues aux conflits sont aussi minimisees par Implication d'un facteur 
beta dans Toperation de marquage. 

Dans un mode de realisation particulier de 1' invention, les forces de 
marquage utilisees sont les suivantes : F R =0,4 (pour la composante rouge) ; 
F G =0,4 (pour la composante verte) ; F B =0,8 (pour la composante bleue), avec 
10 beta = 0,9 (en cas de conflit). 

3.6 Reconstruction de Vintage 

L' image est ensuite reconstruite par transformee en ondelettes inverse, en 
tenant compte des coefficients modifies par le marquage. Ainsi, la figure 8b 
illustre l'image « maison » marquee, reconstruite par transformee en ondelettes 
1 5 inverse, et la figure 8a presente la merae image avant marquage. 
4. Detection de la marque 

On decrit desormais les differentes operations mises en ceuvre lorsque l'on 
cherche a detecter une marque dans une image couleur tatouee selon l'invention. 
En effet, apres traitement de la marque (transmission, compression JPEG, filtrage 
20 median ou ajout de bruit), la marque peut etre detectee. 
4.1 Decomposition en ondelettes 

L'operation de decomposition dans le domaine des ondelettes de l'image 
marquee est identique a celle effectuee dans le cadre du marquage, presentee au 
paragraphe 1 ci-dessus. 
25 4.2 Definition des vecteurs et choix de la convention 

De m§me, les operations de definition des vecteurs et du choix de la 
convention de marquage sont identiques a celles definies precedemment dans les 
paragraphes 3.1, 3.2 et 3.3: 

4.3 Recuperation d'un bit de la marque detectee 
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Pour chaque coordonnee, on verifie la proximite entre le vecteur marque et 
les deux vecteurs de reference. Suivant la convention choisie, il est possible de 
reconnaitre la valeur du bit integre. Le tableau 1 ci-dessous presente tous les cas 
de figure possibles dans la reconstruction de la marque. La marque est ainsi 



5 reconstruite bit par bit. 





convention 1 


convention 2 


convention 3 






< 


v.M,y)-v c {x,y)\ 


0 


0 


1 






< 


PAx,y)-y B {x,y)\ 


1 


1 


0 




V B {x,y)-V A {x,y) 


< 


y a (x,y)-V c (x,y)\ 


0 


1 


1 




V B (x,y)-V c (x,y) 


< 


V B (x,y)-V A (x,y)\ 


1 


0 


0 




V c {x,y)-V A {x,y) 


< 


V c {x,y)-V B {x,y)\ 


0 


1 


0 




V c {x,y)-V B (x,y)\ 


< 


V c {x,y)-V A {x,y)\ 


1 


0 


1 



4 A Reconstruction de la signature 



La signature detectee est ensuite reconstruite a partir de la marque 
detectee. Pour cela, la moyenne des bits de la marque correspondant a une 
coordonnee de la signature est effectuee, selon un principe similaire a celui utilise 
10 lors de la construction de la marque decrite au paragraphe 2. Ainsi, la signature : 
detectee S est reconstruite, bit par bit. 

En d'autres termes, on commence, pour chaque coordonnee des images de 
details, par detecter chaque bit de la marque, en s'appuyant sur la tableau ci- 
dessus, qui indique la valeur du^bit insere, pour chacune des conventions des 
15 figures 6a a 6c, en fonction des distances entre vecteurs. 

Dans un mode de realisation particulier de l'invention, la signature S de 
16 bits est ensuite reconstruite par le calcul de la moyenne des 16 valeurs de la 
marque pour chaque coordonnee correspondante, suivant la forme de la 
redondance definie dans le paragraphe 2. 
20 4.5 Reconnaissance de la signature marquee par rapport a la 

signature initiale 
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II reste a definir le taux de correspondance entre la signature initiale et la 
signature detectee. Pour cela, on effectue le calcul de la correlation entre les deux 
signatures S et S par l'operation suivante : 

5 Les deux signatures correspondent si cc(S,S)^T , ou T est le seuil de 

decision pour la correspondance des signatures. 

Ce calcul de correlation permet de definir le taux de ressemblance entre la 
signature originale S et la signature detectee S . Dans un mode de realisation 
particulier de 1'invention, on fixe le seuil T a 0,7, de sorte que Ton considere que 
10 si cc(S,S) ^ 0,7 , la signature detectee correspond a la signature originale. 

Comme indique precedemment, les figures 8a et 8b presentent une image 
originale, et V image marquee correspondante, selon les parametres du mode de 
realisation particulier de 1'invention expose ci-dessus. 

Les figures 9a a 9c illustrent, sous forme de graphes, les resultats de tests 

15 de robustesse de la marque inseree sur l'image « maison » des figures 8a et 8b. 

Les graphes des figures 9a a 9c correspondent respectivement aux tests de la 
marque face aux attaques de compression JPEG, a l'ajout de bruit et au filtrage 
median, et represented la moyenne de 20 mesures (utilisation de 20 clefs K 
differentes) pour chaque attaque. Comme on peut le voir, les resultats sont stables, 

20 et la robustesse de la technique de tatouage de 1'invention contre ces attaques est 
tres bonne. 
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ANNEXE 1 : RAPPEL SUR LA THEORIE DES ONDELETTES 

1. Introduction 

La theorie des ondelettes permet d'approximer une fonction (courbe, surface, etc.) 
a differentes resolutions. Ainsi, cette theorie permet de decrire une fonction sous 
5 la forme d'une approximation grossiere et d'une serie de details permettant de 
reconstruire parfaitement la fonction originale. 

Une telle representation multi-resolution (Mallat S., "A Theory for Multiresolution 
Signal Decomposition: the Wavelet Representation", IEEE Transactions on 
Pattern Analysis and Machine Intelligence, Vol. 11, No. 7, Juillet 1989, pp. 674- 

10 693.) d'une fonction permet done d'interpreter de fa?on hierarchique 
Information contenue dans la fonction. Pour ce faire, cette information est 
reorganisee en un ensemble de details apparaissant a differentes resolutions. 
Partant d'une sequence de resolutions croissantes (rj) jez> les details d'une 
fonction a la resolution r y sont definis comme la difference d'information entre 

15 son approximation a la resolution r } et son approximation a la resolution r y+1 . 

2. Notations 

Avant de presenter plus en detail les fondements de l'analyse multi-resolution, on; 
presente dans cette section les notations utilisees. - 

❖ Les ensembles des nombres entiers et reels sont respectivement notes Z et 

20 /?. 

♦> L 2 (R)d6note l'espace vectoriel des fonctions unidimensionnelles f(x) 

mesurables et integrables. 

❖ Pour f(x)<EL 2 (R)et g(x)<EL 2 (R), le produit scalaire de/fxj et g(x) est 

defini par : 

25 (f(x),g(x))=jj(u)g(u)du. 

❖ Pour f(x)<EL 2 (R)et g(x)<EL 2 (R), la convolution de f(x) et g(x) est 
definie par : 

/ *g(x)= " )g( x ~ u ) du ■ 
*> L 2 (R 2 ) denote Tespace vectoriel des fonctions f(x, y) de deux variables 
30 mesurables et integrables. 
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❖ Pour f(x,y)<EL 2 (R 2 )ei g(x,y)^L 2 (R 2 ), le produit scalaire dzf(x,y) et 
g(x,y) est defini par : 



Dans cette section, on presente de maniere intuitive les proprietes desirees de 

1'operateur permettant V analyse multi-resolution d'une fonction. Ces proprietes 

sont issues des travaux de Mallat S., dans "A Theory for Multire solution Signal 

Decomposition: the Wavelet Representation", IEEE Transactions on Pattern 

Analysis and Machine Intelligence, Vol. 1 1, No. 7, Juillet 1989, pp. 674-693. 
Soit A 2 j 1'operateur qui approxime une fonction f (x)€zL 2 (R)a la resolution 

2 J (j s> 0) (i.e. f(x) est definie par 2 J echantillons). 
Les proprietes attendues de A oJ sont les suivantes : 

1. A f est un operateur lineaire. Si ^./(x^represente l'approximation de 
f(x) a la resolution 2 j , alors A ;f(x)ne doit pas etre modifie lorsqu'on 

l'approxime de nouveau a la resolution 2 J . Ce principe s'ecrit 
A^ o A 2 j = A v et montre que 1'operateur A 2i est un operateur de projection 

dans un espace vectoriel V 2i <ZL 2 (R). Cet espace vectoriel peut etre 

interprets comrne 1'ensemble de toutes les approximations possibles a la 
resolution 2 J des fonctions de L 2 (R) . 

2. Parmi toutes les approximations possibles de f(x) a la resolution V , 
A 2 jf(x) est la plus similaire a f(x). L'operateur A 2i est done une 

projection orthogonale sur V j . 

3. L' approximation d'une fonction a la resolution 2 J+{ contient toute 
1' information necessaire pour calculer la meme fonction a la resolution 
inferieure 2 J . Cette propriete de causalite induit la relation suivante : 



4. L'operation d'approximation est la meme a toutes les resolutions. Les 
espaces des fonctions approximees peuvent etre derives les uns des autres 
par un changement d'echelle correspondant a la difference de resolution : 




3. 



Proprietes de 1 'analyse multi-resolution 



Vy€EZ,K 2> CF 2 



Vj€iZ,f(x)<EV 2j o /f2ij£K 2 ,. 
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5. Lorsque Ton calcule une approximation de f(x) a la resolution V , une 
partie de Information contenue dans/fx) est perdue. Cependant, lorsque 
la resolution tend vers l'infini, la fonction approximee doit converger vers 
la fonction f(x) originale. De la meme facon, lorsque la resolution tend 
5 vers zero, la fonction approximee contient moins d' information et doit 

converger vers zero. 

Tout espace vectoriel (V 2 , j^qui satisfait I'ensemble de ces proprietes est appele 
approximation multi-resolution de L 2 (R). 

d Analyse multi-resolution d'une fonction unidimensionn elle 

10 4.1 Recherche d'une base de V 2 , 

Comme indique ci-dessus, l'operateur d'approximation ^ est une projection 
orthogonale sur l'espace vectoriel V v . Afin de caracteriser numeriquement cet 
operateur, il faut trouver une base orthonomale de V 2 , . 

K,,etant un espace vectoriel contenant les approximations de fonctions de L 2 (R) 
15 a la resolution 2 y , toute fonction f(x)<EV 2i peut etre vue comme un vecteur a 
2 J composantes. II faut done trouver 2 J fonctions de base. 

Un des theoremes principaux de la theorie des ondelettes stipule qu'il existe une 
fonction unique <$>(x)<EL 2 (R) , appelee fonction echelle, a partir de laquelle on 
peut definir V fonctions de base 9{(x)d& V 2 ,pax dilatation et translation de 

20 <&(x) : 

<T>{ (x) = <$>(2' x-i),i = 0, -- ,2 J - 1 • 
Approximer une fonction f(x)<=L 2 (R) a la resolution V revient done a projeter 
de fa ? on orthogonale f(x) sur les V fonctions de base <$>[(x). Cette operation 
consiste a calculer le produit scalaire def(x) avec chacune des 2' fonctions de 
25 base $>\(x) : 

A 2 J(x) = k ytf(u)M(«))*U x ) 



: \f(u )M2 J u - k))®(2 J u - k). 
On peut montrer que A 2 J(x) peut etre reduit la convolution def(x) avec le filtre 
passe-bas 0>(x) , evaluee au point k : 



1 er depot 



31 



A v f-( f(u) * ^r-2 y u))(k),kGZ. 
Comme <i>(x ) est un filtre passe-bas, A v f peut etre interprets comme un filtrage 

passe-bas suivi d'un sous-echantillonnage uniforme. 
42 Construction de V analyse multi-resolution 
5 Dans la pratique, les fonctions/a approximer (signal, image, etc.) sont discretes. 
Supposons que la fonction f(x) originale soit definie sur n = 2 k (kGZ) 

echantillons. La resolution maximale def(x) est alors n. 

Soit A n fY approximation discrete de f(x) a la resolution n. La propriete de 
causalite (voir §3 ci-dessus) pretend que Ton peut calculer A y f a partir de A n f 

10 pour tout j<k. 

En effet, en calculant la projection des 2 y fonctions de base <t>{(x)de V 2 ,sur 
V Jtl , on peut montrer que A J} f peut etre obtenue en convoluant ^,/avec le 
filtre passe-bas correspondant a la fonction echelle et en sous-echantillonnant le 
resultat d'un facteur 2 : 

15 A 1 J(u) = )i(k - 2u)A 2J .f(k) f 0 <; u < V -1 

avec h(n) ^{^(2u)^(u-n))ynE:Z , 
4.3 La fonction detail 

Comme mentionne dans la propriete (5) du paragraphe 3, l'operation qui consiste 

a approximer une fonction f(x) a une resolution 2 J a partir d'une approximation a 

20 la resolution 2 J+{ genere une perte d' information. Cette perte d'information est 

contenue dans une fonction appelee fonction detail a la resolution 2 J ct notee 
D 2 j f. II est a noter que la connaissance de D^f et A %J f permet de reconstruire 

parfaitement la fonction approximee A^f . 

La fonction detail a la resolution 2 j est obtenue en projetant de maniere 
25 orthogonale la fonction originale f(x) sur le complement orthogonal de dans 

V f4t . Soit W 2 j cet espace vectoriel. 

Pour calculer numeriquement cette projection, il faut trouver une base 
orthonormale de W v ; e'est-a-dire 2 y fonctions de base. Un autre theoreme 

important de la theorie des ondelettes stipule qu'a partir d'une fonction echelle 
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&(x), il est possible de definir V fonctions de base de W v , . Ces fonctions de 
base^Yx; sont obtenues par dilatation et translation d'une fonction Vf(x) 

appelee fonction ondelette : 

W/(x ) = V(2 J x - i), i. = 0,- • -,2 J - 1. 
5 De la meme facon que pour la construction de l'approximation A v f , on peut 
montrer que D 2> /peut etre obtenu par une convolution de la fonction originale 
f{x) avec le filtre passe-haut W(x) suivi d'un sous-echantillonnage d'un facteur 
2 J : 

D 2i f = (f(u)*V>(-2 i u))(k),k<EZ. 

10 4.4 Extension a V analyse multi-resolution de fonction bi-dimensionnelles 

Dans cette section, on presente la maniere d'etendre l'analyse multi-resolution par 
ondelettes aux fonctions de L 2 (R 2 ) comme les images. 

Pour ce faire, on utilise les meme theoremes que ceux utilises precedemment, 
Ainsi, si l'on note F 2 ,l'espace vectoriel des approximations de L 2 (R 2 )a la 

15 resolution 2 y , on peut montrer que Ton peut trouver une base orthonormale de 
V t cn dilatant et translatant une fonction echelle <&(x,y)<EL 2 (R ) : 

&.(x,y) = ®(2 j x - i,Vy - j), (i,j) 6Z 2 . 
Dans le cas particulier des approximations separables de L 2 (R 2 ), on a 
<P(x,y) = ®(x)®(y) oti <p( x ) est U ne fonction echelle de L 2 (R). Dans ce cas, 

20 l'analyse multi-r6solution d'une fonction de L 2 (R 2 ) est effectue en traitant de 
facon sequentielle et separable chacune des dimensions x et y. 
Comme dans le cas unidimensionnel, la fonction detail a la resolution 2 y e st 
obtenue par une projection orthogonale de f(x,y) sur le complement de V v dans 
V 2 j., , note W 2 j . Dans le cas bidimensionnel, on peut montrer que si l'on note 

25 *V(x )\z fonction ondelette associe a la fonction echelle ®(x) , alors les trois 
fonctions definies par : 

W x (x,y) = <$>(x)W(y) 
y 2 (x,y) = y(x)®(y) 

y 3 (x,y) = W(x)Wy) 
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sont des fonctions ondelettes de L 2 (R 2 ) ■ En dilatant et translatant ces trois 
fonctions ondelettes, on obtient une base orthonormale de W 2i : 

l?)(x.y)-<P9(2 J x-k,2 J y-l) 

\y 2 (x,y ) = W(2 J x - k,Vy-l) 

xy>( x ,y ) = W(2 J x - k,2 j y - 1). 

La projection de f(x,y) sur ces trois fonctions de la base de W v donne trois 
fonctions details : 

D\, f = f(x,y) * & (-xWj (-y) 
D 2 2j f = f(x,y)^V j (-x)^/-y) 
D l 2j f = f(x,y)*y J (-x)Wj(-y) 
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REVENDICATIONS 

1. Procede de tatouage d'une image en couleurs, presentant trois 
composantes, 

caracterise en ce qu'il comprend une etape d'insertion d'une marque de tatouage, 
en au moins un point de l'image, selon une regie d'insertion prenant en compte la 
position relative de trois vecteurs de composante, pour chacune desdites 
composantes respectivement, associes audit au moins un point. 

2. Procede de tatouage selon la revendication 1, caracterise en ce qu'on 
selectionne, pour chacun desdits points concernes, deux vecteurs en tant que 
vecteurs de reference et un vecteur a marquer pour porter ladite marque de 
tatouage. 

3. Procede de tatouage selon la revendication 2, caracterise en ce que ledit 
vecteur a marquer est le vecteur central, se trouvant entre lesdits vecteurs de- 
reference. 

4. Procede de tatouage selon la revendication 3, caracterise en ce qu'on 
calcule les distances deux a deux entre lesdits trois vecteurs, lesdits vecteurs de 
reference etant ceux les plus distants l'un de l'autre. 

5. Procede de tatouage selon la revendication 4, caracterise en ce qu'on. 
determine une frontiere entre lesdits vecteurs de reference, definissant deux zones 
associees respectivement aux valeurs binaires « 0 » et « 1 ». 

6. Procede de tatouage selon la revendication 5, caracterise en ce que ladite 
frontiere est definie par la bissectrice entre lesdits vecteurs de reference. 

7. Procede de tatouage selon l'une quelconque des revendications 5 et 6, 
caracterise en ce que le marquage dudit vecteur a marquer comprend un 
deplacement eventuel dudit vecteur a marquer dans l'une desdites zones, selon la 
valeur binaire de marquage a appliquer. 

8. Procede de tatouage selon la revendication 7, caracterise en ce que ledit 
deplacement a une amplitude variable, fonction d'une force de marquage choisie 
en fonction d'au moins une caracteristique locale de ladite image. 

9. Procede de tatouage selon l'une quelconque des revendications 5 a 8, 
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caracterise en ce que le marquage dudit vecteur a marquer comprend un 
deplacement eventuel dudit vecteur a marquer vers Tun ou l'autre desdits vecteurs 
de reference, selon la valeur binaire de marquage a appliquer. 
10. Procede de tatouage selon Tune quelconque des revendications 1 a 9, 

5 caracterise en ce qu'il comprend egalement une etape de transformation en 
ondelettes de chacune desdites composantes de 1' image, et en ce que lesdits trois 
vecteurs de composantes sont determines, pour chaque point d'au moins un 
niveau de decomposition de ladite transformation en ondelettes, pour chacune 
desdites composantes respectivement. 

10 11- Procede de tatouage selon Tune quelconque des revendications 1 a 10, 
caracterise en ce qu'on prevoit au moins deux conventions de marquage d'un 
vecteur. 

12. Procede de tatouage selon la revendication 11, caracterise en ce qu'on 
choisit comme convention de marquage pour une image donnee celle limitant les 

15 risques de conflit lors de la detection de ladite marque. 

13. Procede de tatouage selon la revendication 12, caracterise en ce qu'on 
choisit une desdites conventions en fonction des nombres de vecteurs de reference 
present dans ladite image pour chacune desdites composantes. 

14. Procede de tatouage selon Tune quelconque des revendications 6 a 13, 
20 caracterise en ce que ledit marquage est, au moins dans des situations normales, 

calcule selon Tequation suivante : 

V M , w (^y) = Vr (*>y) - 0 - F M XV, (*,y) - V M {x y y)) 
ou V M (x,y) est ledit vecteur a marquer, V R (x,y) est un desdits vecteurs de 
reference, F M est ladite force de marquage et V M W {x 7 y) est ledit vecteur marque. 
25 15. Procede de tatouage selon les revendications 11 et 14, caracterise en ce 
que, en cas de conflit, ledit marquage est calcule selon l'equation suivante : 

ou beta<L 

16. Procede de tatouage selon Tune quelconque des revendications 10 a 15, 
30 caracterise en ce qu'il comprend, apres ladite etape dissociation d'une marque de 
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tatouage, une etape de transformation en ondelettes inverse, delivrant une image 
marquee. 

17. Procede de tatouage selon Tune quelconque des revendications 1 a 15, 
caracterise en ce que ladite marque est une signature binaire pseudo-aleatoire 
5 ecrite de fagon redondante. 

18* Procede de tatouage selon Tune quelconque des revendications 1 a 17, 
caracterise en ce que lesdites composantes appartiennent au groupe comprenant : 
les composantes RGB ; 
les composantes YUV ; 
10 - les composantes CMY. 

19. Dispositif de tatouage d'une image en couleurs, presentant trois 

composantes, 

caracterise en ce qu'il comprend des moyens d'insertion d'une marque de, 
tatouage, en au moins un point de 1* image, selon une regie d'insertion prenant en^ 
15 compte la position relative de trois vecteurs de composante, pour chacune desdites 
composantes respectivement, associes audit au moins un point. 

20. Programme d' ordinateur comprenant des instructions de code de s 
programme enregistre sur un support utilisable dans un ordinateur pour le; 
tatouage d'une image en couleurs, presentant trois composantes, caracteris^ en ce 

20 que ledit programme comprend des moyens de programmation lisibles par 
ordinateur pour effectuer une etape d'insertion d'une marque de tatouage, en au 
moins un point de l'image, selon une regie d'insertion prenant en compte la 
position relative de trois vecteurs de composante, pour chacune desdites 
composantes respectivement, associes audit au moins un point. 

25 21. Image en couleurs presentant trois composantes et comprenant un tatouage 
obtenu par insertion d'une marque de tatouage, en au moins un point de l'image, 
selon une regie d'insertion prenant en compte la position relative de trois vecteurs 
de composante, pour chacune desdites composantes respectivement, associes audit 
au moins un point. 

30 22. Procede de detection d'un tatouage dans une image tatouee, realise selon le 
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procede de Tune quelconque des revendications 1 a 18, caracterise en ce qu'il 
comprend une etape de recuperation d'une marque de tatouage, en au moins un 
point de Fimage, selon une regie de recuperation prenant en compte la position 
relative de trois vecteurs de composante, pour chacune desdites composantes 
5 respectivement, associes audit au moins un point. 

23. Procede de detection selon la revendication 22, caracterise en ce que ladite 
etape de recuperation comprend une sous-etape de calcul des distances deux a 
deux entre lesdits trois vecteurs, et en ce que les deux vecteurs les plus distants 
Tun de T autre sont des vecteurs de reference, le troisieme vecteur etant le vecteur 

10 marque portant ladite marque de tatouage. 

24. Procede de detection selon la revendication 23, caracterise en ce que, une 
frontiere entre lesdits vecteurs de reference definissant deux zones associees 
respectivement aux valeurs binaires « 0 » et « 1 », ladite etape de recuperation 
comprend egalement une sous-etape de reperage de la zone dans laquelle se situe 

15 ledit vecteur marque, et une etape dissociation d'une valeur binaire 
correspondante. 

25. Procede de detection selon Tune quelconque des revendications 22 a 24, 
caracterise en ce que ladite marque est recuperee au moins deux fois, et en ce 
qu'on met en ceuvre un calcul de correlation par rapport a une signature de 

20 reference, af in de decider si le tatouage est ou non correctement detecte. 

26. Dispositif de detection d'un tatouage dans une image tatouee, realise selon 
le procede de Tune quelconque des revendications 1 a 18, caracterise en ce qu'il 
comprend des moyens de recuperation d'une marque de tatouage, en au moins un 
point de 1'image, selon une regie de recuperation prenant en compte la position 

25 relative de trois vecteurs de composante, pour chacune desdites composantes 
respectivement, associes audit au moins un point. 

27. Programme d'ordinateur comprenant des instructions de code de 
programme enregistre sur un support utilisable dans un ordinateur pour la 
detection d'un tatouage dans une image tatouee, realise selon le procede de Tune 

30 quelconque des revendications 1 a 18, caracterise en ce que ledit programme 
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comprend des moyens de programmation lisibles par ordinateur pour effectuer 
une etape de recuperation d'une marque de tatouage, en au moins un point de 
l'image, selon une regie de recuperation prenant en compte la position relative de 
trois vecteurs de composante, pour chacune desdites composantes respectivement, 
associes audit au moins un point. 
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